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We point out that the third moments of conserved charges, the baryon and electric charge numbers, 
and energy, as well as their mixed moments, change their signs around the QCD phase boundary in 
the temperature and baryon chemical potential plane. These signs can be measured in relativistic 
heavy ion collisions, and will give clear information on the phase structure of QCD and the 
state of the system in the early stage of relativistic heavy ion collisions. The behaviors of these 
moments on the temperature axis and at small quark chemical potential can be analyzed in lattice 
QCD simulations. We emphasize that the third moments obtained on the lattice, together with the 
experimental results, will provide a deep understanding about the QCD phase diagram and the 
location of the state created in heavy ion collisions. 
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1. Introduction 

It is believed that the phase diagram of Quantum chromodynamics (QCD) in the temperature 
(T) and baryon chemical potential (jUb) plane has a critical point where the first order phase tran- 
sition at low T and high ;Ub terminates It is one of the most challenging subject to confirm 
the existence of the QCD critical point experimentally using the relativisitc heavy ion collisions. 
A method suggested for this purpose is to exploit fluctuation observables. The singularity at the 
critical point, at which the transition is of second order, can cause enhancements of fluctuations if 
fireballs created by heavy ion collisions pass near the critical point during the time evolution ^J. 
If such enhancements are sufficiently large and if the fluctuations can survive until the freezeout, 
enhancements of fluctuations will be observed by event-by-event analysis in heavy ion collisions. 
Because of finite size effects and critical slowing down, however, such singularities are blurred and 
its experimental confirmation would be quite difficult BJ]. In fact, so far no clear evidence for 
the critical point has been detected [0]. 

Other proposed way to exploit fluctuation observables to analyze the phase structure of QCD 
are those to use fluctuations of conserved charges |7], ||]. Since some fluctuations of conserved 
charges behave differently between the hadronic and quark-gluon phases, these fluctuations may 
be used as an indicator of the realization of the phase transition. Approaches to use higher or- 
der moments for this purpose have been also suggested recently ^ and experimental attempts to 
measure those higher order moments were reported, for example, in Ref. [|lC|]. 



2. Third moments of conserved charges 

Almost all previous studies focus on the absolute value, especially the enhancement, of each 
observable around the phase boundary. In the present talk, we propose to employ signs of third 
moments of conserved charges around the averages, which we call, for simplicity, the third mo- 



ments in the following, to infer the states created by heavy ion collisions [lip. In particular, we 
consider third moments of conserved quantities, the net baryon and electric charge numbers, and 
the energy, 

m3 , CCC ) 3 «H« <2g>, (2,) 

where N c with c = B,Q represent the net baryon and electric charge numbers in a subvolume V, 
respectively, E denotes the total energy in V, 8N C = N c — (N c ), and 8E = E — (E). We also make 
use of the mixed moments defined as follows: 

m3 ( ccE ),«W>, ^ EE)s <fW>. ( , 2) 

To understand the behaviors of these moments around the QCD phase boundary, we first notice 



that the moments Eqs. (2.1) and (2.2) are related to derivatives of the thermodynamic potential per 
unit volume, ft), up to third order with respect to the corresponding chemical potentials and T. The 
simplest example is ni3(BBB), which is given by 

m 3 (BBB) = -— -j = - — , (2.3) 
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where the baryon number susceptibility, Xb » is defined as 

& = _^ = « (24 ) 

The baryon number susceptibility Xb diverges at the critical point and has a peak structure around 
there [§, [l2, 13]. Since m3(BBB) is given by the /Ib derivative of Xb as in Eq. (2.3), the existence of 



the peak in Xb means that ra3(BBB) changes its sign there. Although the precise size and shape of 
the critical region are not known, various models predict that the peak structure of Xb well survives 
far along the crossover line [|l2|, [[]] (See, the left panel of Fig. [I] as a demonstration of this feature in 
a simple effective model; the details will be explained later). This means that the near (hadron) and 
far (quark-gluon) sides of the QCD phase boundary can be distinguished by the sign of /M3(BBB) 
over a rather wide range around the critical point. As we shall see later, all third moments presented 



in Eqs. (2.1) and (2.2) can be expressed in terms of derivatives of corresponding susceptibilities, 
which diverge at the QCD critical point and hence change their signs there. 

The third moments can be measured in heavy ion collisions by the event-by-event analysis 
similarly to fluctuations, provided that N c and/or E in a given rapidity range, Ay, in fireballs cre- 
ated by collisions is determined in each event. The measurement of Nb is difficult because of the 
difficulty in identifying neutrons. On the other hand, ,/Vq and E can be measured with the existing 
experimental techniques. Four out of the seven third moments in Eqs. (2T) and (^2|) composed of 



A^q and E thus can be determined experimentally. 

All quantities we are considering here, A^Q and E, are conserved charges and the variation of 
their local densities requires diffusion. In Ref. [j7[], it was shown that the effect of diffusion is small 
enough for the fluctuations of the baryon and electric charges if the rapidity range is taken to be 
Ay > 1. In the estimate, the one dimensional Bjorken expansion and straight particle trajectories 
were assumed. If the contraction of hadron phase due to the transverse expansion and the short 
mean free paths are taken into account, the above estimate will be more relaxed. 

Once the negativeness of third moments is established experimentally, it is direct evidence of 
two facts: (1) the existence of a peak structure of corresponding susceptibility in the phase diagram 
of QCD, and (2) the realization of hot matter beyond the peak, i.e. the quark-gluon plasma, in 
heavy ion collisions. We emphasize that this statement using the signs of third moments does not 
depend on any specific models. The experimental measurements of signs of moments also have 
an advantage compared to their absolute values: it is usually essential to normalize experimentally 
obtained values by extensive observables, such as the total charged particle number A^h, in order 
to compare the experimental results with theoretical predictions [0, §]. In the measurement of 
signs, however, normalization is not necessary. It is this feature that our proposal is less subject to 
experimental and theoretical ambiguities and more robust than previously proposed ones. 



3. Analogy to mountain climbing 

The advantage of using signs of third moments, instead of enhancements of absolute values, 
would be nicely explained in an analogy to mountain climbing. Now, we theoretically expect the 
existence of "mountains" of susceptibilities (See, the left panel of Fig. [[]) in the QCD phase dia- 
gram. Experimentalists try to confirm the existence of the peak in heavy ion collisions, in other 
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words, by directly climbing the mountains. During this expedition, in order to confirm the exis- 
tence of the peak previous strategy using fluctuation observables tried to measure the altitude, i.e. 
the value of susceptibility itself. On the other hand, our proposal to use the third moments would 
be compared to the measurement of the derivative of trails. If one has experiences of climbing 
mountains, one knows that the measurements of altitude is difficult, but it is quite easy to recognize 
whether one is going up or down, namely to recognize just the sign of the derivative of the slope, at 
each moment of mountain climbing. In particular, if one has ever arrived at the edge of the moun- 
tains, one knows that it is the most impressive moment during the mountain climbing; derivative 
of the trail changes positive to negative there, and completely different scenary suddenly manifests 
itself in front of you. Such a moment is so unique that one can clearly realize that he/she arrives 
at the edge. To confirm the existence of the moutanins, therefore, the altimeter is not necessary. 
Measuring the derivative of trail is a much easy and robust way for this purpose. 

In terms of the third moments of conserved chages, experimental confirmation of the sign 
should be easier than that of the absolute value. Furthermore, once the negativeness of third mo- 
ments is established experimentally, it is quite strong evidence of the existence of peak structure 
of corresponding susceptibility; no model dependences enter this statement. Of course, measuring 
the negative third moments in experiments depends on whether the fireballs can remember the im- 
pressive moment until the freezeout, or not. Exploiting the moments of conserved chages plays a 
crucial role for this discussion, as argued above. As the experience of arriving at the edge would be 
so impressive, why don't we anticipate that the we can remember the moment even after you come 
back home? 



4. Other third moments 

Let us now consider the behavior of third moments other than /M3(BBB) around the critical 
point. First, the third moment of the net electric charge wi3 (QQQ) is calculated to be 



d 3 co ld 3 (0 3 d 3 co 3 d 3 co \d 3 co 
m 3 QQQ =-—3= -7TT-«IT-»3 2-qITT' ( 4 -!> 



where /Iq represents the chemical potential associated with Nq, i.e. d/d/j.Q = (2/3)d/djU u — 
(l/3)<?/djUd = (d/d^lB + d / dpL\) /2, and the isospin chemical potential is defined as pL\ = (/x u — 
/id) /2 with ju u ,d being the chemical potentials of the up and down quarks, respectively. In relativis- 
tic heavy ion collisions, the effect of isospin symmetry breaking is small. Assuming the isospin 



symmetry, the second and last terms in Eq. (4.1) vanish and one obtains 



i»3(QQQ) = ^0fc + 3^) s (4.2) 

with the isospin susceptibility %\ = —d 2 (o/dnf. Under the isospin symmetry, %\ does not diverge 
at the critical point because the critical fluctuation does not couple to the isospin density [Bj. The 



critical behavior of the term in the parenthesis in Eq. (4.2) in the vicinity of the critical point is 
thus solely governed by Xb- Since niT,(QQQ) is a /Ib derivative of this term, a similar behavior as 
ra3(BBB) is expected. 



4 



Third moments of conserved charges 



Masakiyo Kitazawa 



Next, it can be shown that mixed moments including a single E are concisely given by 

1 d(T X c 



mi, (ccE) 



(4.3) 



dT 

with c = B, Q, where %q = —d co/dHq = {%b + Xi)/4 is the electric charge susceptibility. The T 
derivative in Eq. ( fO| ) is taken along the radial direction from the origin with fixed jX = Hb/T, i.e. 



d/dT\fi = d/dT\n B + (fj,^/T)d /d/ielr- Since T% c diverges at the critical point, Eq. (4.3) again 
leads to a similar behavior of m^^ccE) as the above-mentioned moments. 

To argue the behaviors of remaining third moments including two or three E's, it is convenient 
to first define Cfx = —T{d 2 co/dT 2 )fx = {(8E) 2 ) /VT 2 . The third moments are then given by 



1113 (EEE) 



i <?(r 2 c A ) 



r 3 dT 



1 c)C ' 

m 3 (BEE) = 2m 3 (QEE) = -— £ . (4.4) 



Since is the second derivative of ft) along a radial direction, it diverges at the critical point 
which belongs to the same universality class as that of the 3D Ising model. Therefore, m3(EEE), 
ra3(BEE), and m 3 (QEE), all change their signs at the critical point. 

5. Region with negative third moments and possible lattice analysis 

While the above arguments, based on the divergence of second derivative of ft), guarantee the 
appearance of the region with negative third moments in the vicinity of the critical point, they do 
not tell us anything about the size of these regions in the T-\Ib plane. In fact, all third moments 
considered here become positive at sufficiently high T and jU B > 0, where the system approaches a 
free quark and gluon system. The regions are thus limited more or less near the critical point. 

The information about the behavior of the third moments at small /Ib can be extracted from 
the numerical results in lattice QCD. For example, with the Taylor expansion method the ther- 
modynamic potential is calculated to be ft) = —C2iT)\i\ — c 4 {T)ji.^ — c^{T)}X^ — • • • , and one can 
read off the behavior of m 3 (BBB) at small /i B as m 3 (BBB) = 24{c 4 (T)n B + 5c 6 (r)/i| H ]. Lat- 
tice simulations indicate that c^{T) is positive definite, while C(,(T) becomes negative in the high 



temperature phase [ 15]. From this result one sees that ra3(BBB) is positive for small /^b, while the 
negative C(,(T) suggests that the sign of »i3(BBB) eventually changes at sufficiently large ;Ub- Other 
moments for small /^b can also be evaluated in the Taylor expansion method by expanding ft) with 
respect to T and ;Uq. If the contour lines of vanishing third moments are close enough to the T-axis, 
the lattice simulations may be able to determine these lines. Since the region with a negative third 
moment should depend on channels, combined information of signs of different third moments, 
and the comparison of the third moments obtained by experiments and lattice simulations, will 
provide a deep understanding about the state of the system in the early stage of relativistic heavy 
ion collisions and the QCD phase diagram. 

6. Analysis in a toy model 

The range of Hb/T where lattice simulations are successfully applied is, however, limited to 
small Hb/T with the present algorithms. In particular, thermodynamics around the critical point 
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Figure 1: Left: T and /Xb dependence of the baryon number susceptibility Xb multiplied by T in the Nambu- 
Jona-Lasinio model. The bold line on the bottom surface shows the first order phase transition line and the 
point at the end is the critical point. Right: Regions where third moments take negative values in the T-jj.Q 
plane. The regions are inside the boundaries given by the lines. 



cannot be analyzed with the Taylor expansion method. In order to evaluate the qualitative behavior 
of the third moments in such a region, one has to resort to effective models of QCD. To make such 
an estimate, here we employ the two-flavor Nambu-Jona-Lasinio model [16] with the standard 
interaction J£ mX = G{(\j/\j/) 2 + (l/Z/ysT,-^) 2 }, where y denotes the quark field. For the isospin 
symmetric matter, this model gives a first order phase transition at large /Ib, as shown on the 
bottom surface of the left panel of Fig. [j] by the bold line. 

In the left panel of Fig. [j], we also show the T and /Ib dependence of T% B calculated in 
the mean-field approximation. One observes that Xb diverges at the critical point, and the peak 
structure well survives along the crossover line up to higher temperatures [p^]. The region where 
each moment becomes negative in the T-}Xb plane is shown in the right panel of Fig. [j]. The figure 
shows that areas with wi3(BBB) < and /M3(BBE) < extend to much lower /Ib and much higher 
T than the critical point. This suggests that even if the critical point is located at high jUb the 
negative third moments can be observed by heavy ion collision experiments. The figure also shows 
that the areas have considerable thicknesses along the radial direction. Since the system stays near 
the phase transition line considerably long regardless of the order of the phase transition, first order 
or crossover, once the state on the far side is created, negative third moments are very likely to 
be formed and observed. The wide regions of negative moments also indicate that they are hardly 
affected by critical slowing down during the dynamical evolution of fireballs. 

The right panel of Fig. [| also shows that areas with negative /713 (EEE) , /713 (QEE) , and mj, (BEE) 
are much larger than those of the other moments in the T-/Ib plane; although not shown in the fig- 
ure, these areas extend even to the T-axis. The behaviors of (EEE) and 7713 (cEE) near the T-axis 
can be checked directly by the lattice simulations. If the range of T satisfying wi3(EEE) < is 
sufficiently wide at /Ib = 0, it is possible that the negative third moments are measured even at the 
RHIC and LHC energies. Whether the negative moments survive or not in this case depends on 
the diffusion time of the energy density, in other words the heat conductivity. One can thus use the 
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signs of m3 (EEE) and 1113 (cEE) to estimate the diffusion time of the charges and energy. The third 
moments m3(QQQ) and m3(QQE), on the other hand, become negative only in a small region near 



the critical point. These behaviors come from the large contribution of %\ in Eq. (|4.2|). 
7. Summary 

In this talk, we have pointed out that the third moments of conserved charges, the net baryon 
and electric charge numbers and the energy, change signs at the phase boundary corresponding to 
the existence of the peaks of susceptibilities. If the negative third moments grow at early stage 
of the time evolution of fireball created in the collisions and if the diffusion of charges is slow 
enough, then the negative third moments will be measured experimentally through event-by-event 
analyses. Once such signals are measured, they serve as direct evidence that the peak structure 
of corresponding susceptibility exists in the phase diagram of QCD, and that the matter on the far 
side of the phase transition, i.e. the quark-gluon plasma is created. The combination of the third 
moments of different channels, and their comparison with the numerical results in lattice QCD 
will reveal various issues on the phase structure and initial states created in heavy ion collisions at 
different energies. 

References 

[1] M. A. Stephanov, PoS LAT2006, 024 (2006). 

[2] M. A. Stephanov, K. Rajagopal, and E. V. Shuryak, Phys. Rev. Lett. 81, 4816 (1998); Phys. Rev. D 60, 
114028(1999). 

[3] Y. Hatta and M. A. Stephanov, Phys. Rev. Lett. 91, 102003 (2003) [Erratum-ibid. 91, 129901 (2003)]. 
[4] B. Berdnikov and K. Rajagopal, Phys. Rev. D 61, 105017 (2000). 
[5] C. Nonaka and M. Asakawa, Phys. Rev. C 71, 044904 (2005). 

[6] I. Arsene et al. [BRAHMS Collaboration], Nucl. Phys. A 757, 1 (2005); B. B. Back et al. [PHOBOS 
Collaboration], Nucl. Phys. A 757, 28 (2005); J. Adams et al. [STAR Collaboration], Nucl. Phys. A 
757, 102 (2005); K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 757, 184 (2005). 

[7] M. Asakawa, U. W. Heinz, and B. Muller, Phys. Rev. Lett. 85, 2072 (2000). 

[8] S. Jeon and V. Koch, Phys. Rev. Lett. 85, 2076 (2000). 

[9] S. Ejiri, F. Karsch, and K. Redlich, Phys. Lett. B 633, 275 (2006); M. A. Stephanov, Phys. Rev. Lett. 
102, 032301 (2009); M. Cheng et al, arXiv:081 1.1006 [hep-lat]. 

[10] T. K. Nayak [STAR Collaboration], arXiv:0907.4542 [nucl-ex]. 

[11] M. Asakawa, S. Ejiri and M. Kitazawa, arXiv:0904.2089 [nucl-th]. 

[12] Y. Hatta and T. Ikeda, Phys. Rev. D 67, 014028 (2003). 

[13] H. Fujii, Phys. Rev. D 67, 094018 (2003). 

[14] C. Nonaka, B. Muller, S. A. Bass, and M. Asakawa, Phys. Rev. C 71, 051901 (2005). 

[15] C. R. Allton et al, Phys. Rev. D 71, 054508 (2005). 

[16] T. Hatsuda and T. Kunihiro, Phys. Rept. 247, 221 (1994). 



7 



